A=eye(4); B=speye(4); [i,j]=find(A);
S1=A(i,j); whos
S1=B(i,j); whos
A(i,j)=2*B(i,j); whos
B(i,j)=A(i,j)^2; whos
